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Abstract 

We consider interaction effects in a granular normal metal at not very low 
temperatures. Assuming that all weak localization effects are suppressed by 
the temperature we replace the initial Hamiltonian by a proper functional 
of phases and study the possibility for a phase transition depending on the 
tunneling conductance g. It is demonstrated for any dimension that, while 
at small g the conductivity decays with temperature exponentially, its tem- 
perature dependence is logarithmic at large g. The formulae obtained are 
compared with an existing experiment and a good agreement is found. 
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In an experimentS on granular Al-Ge thick films several interesting effects have been 
discovered. Destroying superconductivity by a magnetic field up to 17T the authors could 
study, in particular, properties of the normal state. Some features of the normal state related 
to a negative magnetoresistance due to superconducting fluctuations have been discussed 
recentlya but an unusual observation remained completely unexplained. 

What we have in mind is a peculiar temperature dependence of the conductivity found in 
some samples. Samples that had a high room temperature resistivity showed an exponential 
decay of the conductivity as a function of temperature. This behavior is typical for insulators 
and has been interpreted in Ref.0 in this way. Samples with larger intergranular couplings 
did not show any exponential decay but the resistivity did not saturate at low temperatures 
and the authors described its temperature dependence by a power law 

R = AT~ a , (1) 

with a = 0.117. Apparently, with such a small value of a a logarithmic temperature depen- 
dence 

R = A(l-a\nT) (2) 

(obtained by expansion of Eq. ([]]) in a) could describe the experimental data as well. It is 
relevant to emphasize that the array of the grains was three dimensional and one could not 
attribute such a behavior to weak localization effects. 

In this paper, we consider a model for granular metals at not very low temperature and 
demonstrate that changing the dimensionless tunneling conductance g one can have either 
exponential temperature dependence of the resistivity at small g < 1 or the logarithmic 
behavior, Eq. (fj), at large g > 1. It will be shown that the result is applicable for any 
dimensionality of the array of grains, which contrasts usual logarithmic corrections due to 
interference effects! typical for 2D. The Hamiltonian describing the model is chosen as 

H = H + H t + H c , (3) 

where Hq is the one-electron Hamiltonian of isolated grains including disorder within the 
grains. The tunneling of the electrons between the grains is given by 

Ht = tijftada'j- (4) 

i,j,a,a' 

where the summation is performed over the states a, a' of each grain (spin is conserved) 
and over neighboring grains i and j. The possibility to tunnel from the state a to an 
arbitrary state a' of other grains introduces an additional disorder resulting in a finite 
tunnel conductance. 

The term H c in Eq. (|3|) describes the charging energy 

H c = ^J2^C^N r (5) 
Z ij 

In Eq. @, Aj = J2a I $a ( r i) ^Pa ( r i) d r i — N is the excess number of electrons in the z-th 
grain. (N is the dimensionless local potential) and is the capacitance matrix. Eq. ([S]) 
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describes the long range part of the Coulomb interaction in the limit of weak disorder inside 
the grains and has been used in many works. Calculations with the Hamiltonian H, Eqs. 
(||, |), can be replaced by computation of a functional integral over anticommuting fields 
i) a (r). 

Although the model described by Eqs. (|3], |5|) contains only the long range part of the 
Coulomb interaction, it is still very complicated, because at very low temperatures inter- 
ference becomes very important and one has to consider an interplay of localization and 
interaction effects. One could do this either using diagrammatic expansions! or writing 
a non-linear cx-modeli. Both methods allow to consider the limit of large tunneling con- 
ductances g and the results are strongly dependent on the dimensionality. However, the 
behavior, Eq. (fj) or Eq. (||D was not predicted for 3D in any of these works. 

The model, Eqs. ([3], ||), becomes simpler if the temperature T is not very low such 
that low energy diffusion modes are damped. As it was discussed in a recent publication!, 
the granular metal can be described at temperatures T ^> g8, where S is the mean level 
spacing in a single grain, by the Ambegaokar, Eokern and Schon (AES)l functional ot the 
free energy. If g ^ 1, this condition should be replaced by T ^> 5. The limit of not very low 
temperatures not only simplifies the consideration but is interesting on its own because it 
leads to an unusual behavior of physical quantities and is easily accessible experimentally. In 
particular, we will see that changing the tunneling conductance g one may have a transition 
from the exponential temperature dependence of the resistivity to the logarithmic behavior, 
Eq- ©. 

We calculate the conductivity a [uS) using the Kubo formula and making an analytical 
continuation from Matsubara frequencies iu n to real frequencies In order to reduce the 
calculation of physical quantities to a computation of correlation functions with the AES 
action we decouple the interaction term, Eq. (|5|), by integration over an additional Vi (r) 
and then, following Refs.H!, remove this field from Hq by the gauge transformation 

^ (r) - e-^«V< (r) , & (r) = V (r) . (6) 

This is not a trivial procedure, because the new fields ip a (r) should obey, as before, the 
boundary condition ip a (r) = — ip a (r + (3) , (3 = 1/T. So, the field Vi (r) cannot be com- 
pletely removed from Hq at arbitrarily low temperatures. Instead of (pi(r) let us consider 
Mr): 

fa (r) = (f>i (r) + 2nTk i r, (7) 

where — oo < (r) < oo, (pi (0) = (pi and ki = 0, ±1, ±2, .... are the so called winding 
numbers. Performing the gauge transformation with (pi instead of (pi(r), the antiperiodicity 

of the ip a is preserved, but the action still contains a term linear in i(Vi — </>J G (— wrT, mT). 
Only in the limit T ^> 5 this term can be neglected. The integration over 0, (r) implies 
integration over (pi (r) and summation over ki. At large g 3> 1, one can put all ki = 0. 
However, at g ^ 1 one should sum over all ki and neglecting the contribution of the 
non-zero winding numbers, as done in Ref.l, leads to incorrect results. 

Using the phase representation one can write the conductivity o (u) in the form 

,d-2 [ r 1 

(8) 



la 

a{u) = 

u> 



dre tQnT K 



-iuj+8 
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K(r) = (X 2 a (r))-5:(X 1 a (r)X 1 a i (0)), 

i 



X 2 a (r) = e 2 7cg dr' (5 (r) - S (r' - r)) a (r') cos (0 M+a (r') - iii+a (0) 

-^5 ( r ) = e7T 9 a i r ~ t') sin (0i, i+a (r') - 0i,i+ a (r)) , a (r) = T 2 (Re (sin (vrTr + 

where a is a vector connecting the centers of neighboring grains i and i + a, a = |a| , and d 
is the dimensionality of the array. In Eqs. (|8|), 0ij (r) = (pi (r) — (f)- } (r) for i and j standing 
for neighboring grains and 

(...) = J (...) exp (-S) D4> (J exp (-S) D^j ' , (9) 

where Dip stands for both the functional integration over (r) and summation over the 
winding numbers ki. The AES action S can be written as 

S = S c + S t , (10) 

where S c describes the charging energy 

S c = ±Y,^drc /^ T)d ^ T) (11) 
2e 2 Y^o 3 dr dr K J 

and S t stands for tunneling between the grains 

S t = 2ng £ ["drdr'a(r-r') S ^( ^ {T) -^ {T,) ) (12) 



|i-Jl=a 



The dimensionless conductance g is given by g = 2iru 2 t^, where is the tunneling amplitude 
from grain i to grain j (spin is included). 

Although the model described by Eqs. (|B|-|T2]) is simpler than the initial model, Eqs. 



, , explicit formulae can be written only in limiting cases. The same action S, Eqs. (|10|- 
|T2"D , was used in Ref.i, and a metal- insulator transition has been predicted in a 2D array 
of tunnel junctions. However, the authors of Ref.i did not calculate the conductivity but 
discussed properties of the partition function. For large g they did not account for phase 
fluctuations properly which, as we show here, are responsible for the behavior, Eq. (^|). 
Moreover, we find a transition in any dimensionality. 

If the temperature T is very high, T>£ c ~ e 2 C^ , where E c is the electrostatic energy 
of adding one electron to a grain, fluctuations of the phases (j) are negligible and one can set 
(p = in the expressions for Xi and X 2 in Eqs. (§). Then, we obtain easily the conductivity 

do = e 2 ga d - 2 , (13) 

which shows that at such temperatures charging interactions are not important. 
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In the opposite limit T <C E c , transport in the granulated system has much more interest- 
ing characteristics. This inequality can be compatible with the inequality T 3> max{gS, 5}, 
used for the derivation of Eqs. because E c ^$> 5 for 2D and 3D grains. 

We calculate the conductivity at temperatures T <C E c in the limits g ^> 1 and g ^ 1. 

In the limit of large conductances g ^> 1, fluctuations of are small and all non-zero 
winding numbers ki can be neglected. Non-zero ki (as well as variations of iV, ) would lead 
to contributions of order exp (—g) , and can be neglected in any expansion in 1/g. 

Keeping only quadratic in <p terms in Eqs. (|K}|12|) we reduce the action S to the form 



S = T£< n G->_ q ,_ n , (14) 



q,n 



Gqi = co 2 J (4£ (q)) + 2g M £ i 1 ~ cos qa) > 

a 

where (q) = e 2 / (2C (q)) and C (q) is the Fourier-transform of the capacitance matrix CV,- 
(q are quasi-momenta for the array of the grains). One should sum in Eq. (H) over d unit 
lattice vectors a, where d is the dimensionality of the array. 

Keeping only quadratic in <p terms in the action but not expanding the function X 2 , Eq. 
(|), one reduces the correlator (X 2a (^n)) to the form 

(X 2a (Q n )) = 7ie 2 g a (r) (l - e iQ " T ) e~^dr, 



G a (r) =4Ty _ _jG qn siii — sin — (15) 

One can check that the contribution coming from the correlator of the functions X\ in 
Eqs. (§) contains additional powers of 1/g and can be neglected in the main approximation. 
It is very important that the correlator (X\X\) in Eqs. (f|) contains a summation over j, 
which corresponds to the zero quasi-momentum of the function K. If we carried out the 
computation for a single grain the contribution from (X\X\) would not be smaller than the 
one from (X2). 

What remains to be done in order to calculate the conductivity for g ^> 1 is to compute 
the integral in Eq. fli"5] ) for the Matsubara frequencies Q n and make the analytical continu- 
ation Q n — > — iu + 5. In the lowest order in a the result for the conductivity a in the limit 
uj — > is 

a = a (1 - a In (gE c /T)) , a = {^ngdy 1 . (16) 

Thus, at large j > 1 the conductivity decays with temperature logarithmically. Of 
course, one may not use this formula at very low temperature because the present consid- 
eration is valid at not very low temperatures T ^> gS when the AES action may be used. 
At lower temperatures, one should take into account interference effects and, depending on 
the dimensionality d of the array, both metal and insulating states are possible. In contrast, 
Eq. fllBD, is valid in any dimensionality. 
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If we used Eg. (|15D exactly we would obtain the power law, Eq. (|l|), with the exponent a 
from Eq. (|16|) . A similar dependence was written for the voltage dependence of the conduc- 
tance of a single junction in a model with an electromagnetic environment0'0. However, 
taking into account non-quadratic terms in the expansion of the action in phases <fi changes 
this result and one comes again to Eq. ([16]). For a single junction this result is known since 
the works of Refill where a proper renormalization group (RG) equation was written. One 
can check that, in the first order of the RG, the equation is the same for a granular metal 
and does not depend on the arrangement of the grains in the array. It can be written as 

d 9 (0 = L (17) 

d£ 2nd { ' 

where £ = — In r and g (£) is the effective conductance. 

Solving Eq. (|17]) with the boundary condition g (0) = g we come immediately to Eq. 

©• 

However, next orders in the RG (expansion of the Gell-Mann-Low function in g (£)) are 
dependent on the type of the array and differ from the corresponding terms for a single 
junction. The applicability of the one- loop approximation, Eq. (O), and of its solution, 
Eq. (0), implies an additional inequality for the temperature: T > T c = gE c exp(—l/a). 
In other words, Eq. ( |T6"D remains valid until a/e 2 becomes of the order unity. For not very 
large grains this inequality is less restrictive because T c is exponentially small for large g. The 
dimensionality of the array d enters Eqs. (|TB|, |T7|) as a parameter only. Actually, the number 
of contacts with neighboring grains (coordination number) rather than the dimensionality 
itself enters Eqs. ([L6|, |17]). This difference may be important in situations when the grains 
are close packed in a cubic lattice. 

The logarithmic behavior, Eq. (|T6|) , describes the granular system at sufficiently large 
g ^ 1. At smaller g, the temperature dependence becomes exponential and we check this 
statement in the limit g< 1 expanding the functional integral in Eq. (|9|) in the tunneling 
part S t , Eq. flT2D, of the action. The main contribution comes again from the function 
(X 2 (t)) in Eqs. @. In the lowest order one can completely neglect S t , which leads to 
computation of the correlator n (r) 

H (r) = (exp (-< (fa (r) - fa (0))))^ , (18) 

where the phases fa (r) are introduced in Eq. (|7|), S c is given by Eq. ( |TTD and the averaging 
should be performed with this functional (Strictly speaking the function n (r) is sufficient for 
calculating K (r) only for diagonal Cy. However, a proper modification for an arbitrary CV,- 
is simple). The computation of the average in Eq. fll8|) can be performed using two different 
methods. A more straightforward way of calculating is to use the definition of fa (r), Eq. 
(0), which allows to represent the action S c as S c = S c [<j>] + S c [k] and carrying out integration 
over the phase </> and summation over the winding numbers separately. Integrating over the 
phase fa (r) we obtain for < r < (3 (see alsoi!) 



{e -i {MT) - M0))) = ^_ Bii ^ _ Tr 2^ ^ (lg) 

where B l3 = ± (C^V .. 
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However, Eq. fllTf ) is not the final result and the summation of the winding numbers is 
essential. This can be performed using the Poisson summation formula. As a result, the 
function II (r) can be represented as 

n(r) = — e~ B " T e~^ 2TnkBki ~ P ^ k J Bkinkni , (20) 

^ {n k } 

where all are integers and Z is a normalization coefficient (11(0) = 1). The necessary 
periodicity in r of the function II (r) with the period f3 is evident from Eq. (El]) . 

The second method is to use the standard quantum mechanical formalism instead of 
calculating the functional integrals in Eq. fllSD , which has been suggested in an earlier work 
on granulated superconductorslij. Within this approach one writes instead of the action S c , 
Eq. (|TTD , an effective Hamiltonian H e g 

H e s = J2 B ijPiPj, P% = -id/dfa, (21) 

and calculates the thermodynamic average with H e a- For the function IT (r) one should 
calculate the average 

( e -<^)-^°)))£ eff , 0, (r) = e^fre-"^. 

Eigenvalues of the operators pi are integers (eigenfunctions of if e ff must be periodic in with 
the period 2n) and one comes easily to Eq. (|20|) . This consideration explicitly demonstrates 
that accounting for the winding numbers leads to the charge quantization. 

We see from Eq. (^0]) that, in order to get an explicit expression for the conductivity, 
one should sum over all configurations of charge. At high temperatures T ^> E c the sum 
over rifc in Eq. (|20"D can be replaced by integrals and we get II (r) = 1, which leads to Eq. 

(0)- 

In the opposite limit, T -C E c , the main contribution comes from charge configurations 
with the lowest energy. The ground state with all = does not contribute to the 
conductivity. If the lowest excited state corresponds to one charged grain with the charge 
±1 (depending on a particular Cy) we come using Eqs. (||) in the limit u — > to a rather 
simple formula 

a = 2a Q exp(-B u /T). (22) 



Since Ba is the energy corresponding to the charge ±1, Eq. ( P2"D corresponds to conduc- 
tion of an activated electron and hole (the factor 2 means that both of them are taken into 
account). 

Comparing Eq. (p2|) with Eq. ([16]) we come to the conclusion that there must be a 
critical value g c separating in the limit T — > the logarithmic behavior from the exponential 
one. Whether the activation energy (Coulomb gap) turns to zero or has a jump at g — g c 
is not clear from the present consideration. (Strictly speaking, there should not be any 
singularity of the conductivity at finite temperatures but the change of the behavior may be 
noticeable experimentally or numerically). The model of the granular metal may be used 



7 



to describe disordered electron systems at low electron density. In this case, potential wells 
would correspond to the grains. 

The sample of the experimentS that showed the "power law behavior", Eq. (|1|), had 
the room temperature resistivity Ro = 7.3 x 10~ 3 ficm. The diameter of the grains was 
120 ± 20A, which allows, using the value h/e 2 = 4.1 x 10 3 Q, to estimate the dimensionless 
tunnel conductivity as g = 0.7. If we put d = 2 in Eq. fll~6| ) we obtain a = 0.116, which 
exactly corresponds to the experimental value from Eq. (fj). However, the arrays used in 
Ref.0 were rather thick and, at first glance, one should use d = 3. Nevertheless, the value 
of d in Eq. (IB) corresponds rather to the half of the contacts of a single grain than to the 
real dimensionality. Then, the experimental value of a indicates that either the grains are 
not closely packed such that the typical number of contacts per grain is 4 or our calculation 
is too rough to provide a quantitative agreement with the experiment (the value of a, Eq. 
(fL6D, is based on the assumption j > 1 but the experimental value of g is of order 1). 
The resistivity of samples with a high room temperature resistivity behaved as exp (a/T 1 / 2 ^ 
rather than obeying the activation law, Eq. (|22TL But this can be attributed to a variation 
of the size of the grains or of the local potential". 

In conclusion, we suggested a scheme of calculating the conductivity of a granular metal 
at not very low temperatures. On the basis of explicit results we demonstrated the existence 
in any dimensionality of a transition between states with exponential dependence of the 
conductivity on temperature and a logarithmic one. Relating the coefficient a to the room 
temperature conductivity we were able to compare our results with an existing experiment 
and got a good agreement. The model of the granular metal may serve also as a good 
description of disordered systems with a low electron density. 
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